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Abstract 



We summarize the top-quark mass measurements from the CDF and D0 experiments at 
Fermilab. We combine published Run-I (1992-1996) measurements with the most recent 
prehminary Run-II (2001-present) measurements using up to 1 fb~^ of data. Taking 
correlated uncertainties properly into account the resulting preliminary world average 
mass of the top quark is Mt = 170.9 it l.l(stat) it 1.5(syst) GeV/c^, assuming Gaussian 
systematic uncertainties. Adding in quadrature yields a total uncertainty of 1.8 GeV/c^, 
corresponding to a relative precision of 1.1% on the top-quark mass. 



^The Tevatron Electroweak Working Group can be contacted at tev-ewwg@fnal.gov. 
More information can be found at http : / / tevewwg . f nal . gov, 
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1 Introduction 



The experiments CDF and D0, taking data at the Tevatron proton-antiproton coUider located 
at the Fermi National Accelerator Laboratory, have made several direct experimental measure- 
ments of the top-quark pole mass, Mt. The pioneering measurements were based on about 
100 pb"^ of Run-I (1992-1996) data [Il[3[3ll2[5l[6l[71[8l[9l[T0l[IIl[l2]and include results from 
the tt — >• qq'bqq'b (all-j), the tt — > iuqq'bb (l-l-j), and the ti i^ubi'iTb (di-1) decay channel^. 
Results using approximately 350 pb~^ of Run-II (2001-present) data have been published in the 
l-|-j and di-1 channels [T3|, [HI |15l |16]. More recently results using about 1 fb~^ of Run-II data 
have also been published [T71 [18]. The Run-II measurements summarized here are the most 
recent results in the l-|-j, di-1, and all-j channels using 700 — 1000 pb~^ of data and improved 
analysis techniques [HI [HI [191 [201 [211 [22] • 

This note reports the world average top-quark mass obtained by combining five published 
Run-I measurements [21 [3, [SI [3, [IDl HI] with two published Run-II CDF results [171 [IB] , two pre- 
liminary Run-II CDF results [13 [20] and two preliminary Run-II D0 results [211 [22] ■ The com- 
bination takes into account the statistical and systematic uncertainties and their correlations 
using the method of references [231 [21] and supersedes previous combinations [23 [211 [23 [28] . 
The most precise individual measurements of Mt are now the preliminary measurements in the 
1+j channel from Run II. These are 170.9 ± 2.5 GeV/c^ (CDF, [19]) and 170.5 ± 2.7 GeV/c^ 
(D0, [21] )■ These have weights in the new Mt combination of 39% and 40%, respectively. 

The input measurements and error categories used in the combination are detailed in Sec- 
tion [2] and [3], respectively. The correlations used in the combination are discussed in Section [H 
and the resulting world average top-quark mass is given in Section [5] A summary and outlook 
are presented in Section [61 



2 Input Measurements 

For this combination eleven measurements of Mt are used, five published Run-I results, and 
two published plus four preliminary Run-II results. In general, the Run-I measurements all 
have relatively large statistical uncertainties and their systematic uncertainty is dominated by 
the total jet energy scale (JES) uncertainty. In Run-II both CDF and D0 take advantage 
of the larger tt samples available and employ new analysis techniques to reduce both these 
uncertainties. In particular the JES is constrained using an in-situ calibration based on the 
invariant mass of W ^ qq' decays in the l-|-j and all-j channels. The Run-II D0 analysis in 
the l-|-j channel constrains the response of light-quark jets using the in-situ W — >• qq' decays. 

^Here £ = e or /i. Decay channels with exphcit tau lepton identification are presently under study and are 
not yet used for measurements of the top-quark mass. 
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Residual JES uncertainties associated with rj— and pr-dependencies as well as uncertainties 
specific to the response of 6-jets are treated separately. Similarly, the Run-II CDF analysis 
in the 1+j and all-j channels also constrain the JES using the in-situ W —>■ qq' decays. Small 
residual JES uncertainties arising from rj— and pr-dependencies and the modeling of 6-jets are 
included in separate error categories. The Run-II CDF di-1 measurement uses a JES determined 
from external calibration samples. Some parts of the associated uncertainty are correlated with 
the Run-I JES uncertainty as noted below. 

In previous combinations the Run-II CDF 1+j analysis used the JES determined from the 
external calibration as an additional Gaussian constraint. This required us to treat that mea- 
surement as two separate inputs in the combination in order to accurately account for all the 
JES correlations. This Gaussian constraint is not used in the present analysis as it does not 
significantly improve the sensitivity. Thus we can treat this measurement as a single input in 
the same manner as all the other measurements. 

A new analysis technique from CDF is also included (Ixy). This measurement uses the 
mean decay-length from B-tagged jets to determine the top-quark mass. While the statisti- 
cal sensitivity is not nearly as good as the more traditional methods, this technique has the 
advantage that since it uses only tracking information, it is almost entirely independent of 
JES uncertainties. Additionally, since it does not require a full event reconstruction, it can 
use a more inclusive sample of tt candidates (e.g. events with > 3 jets). As the statitistics 
of this sample continue to grow, this method could offer a nice cross-check of the top-quark 
mass that's largely independent of the dominant JES systematic uncertainty which plagues the 
other measurements. The statistical correlation between this measurement and the Run-II CDF 
l-|-j measurement is determined using Monte Carlo signal-plus-background psuedo-experiments 
which correctly account for the sample overlap and is found to be consistent with zero (to 
within < 1%) independent of the assumed top-quark mass. 

The inputs used in the combination are summarized in Table [1] with their uncertainties 
sub-divided into the categories described in the next Section. The correlations between the 
inputs are described in Section HI 



3 Error Categories 

We employ the same error categories as used for the previous world average ^28j. They include 
a detailed breakdown of the various sources of uncertainty and aim to lump together sources of 
systematic uncertainty that share the same or similar origin. For example, the "Signal" category 
discussed below includes the uncertainties from ISR, FSR, and PDF - all of which affect the 
modeling of the ti signal. Additional categories are included in order to accommodate specific 
types of correlations. For example, the jet energy scale (JES) uncertainty is sub-divided into 
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Run-I published 


Run-II preliminary 






CDF 




D0 




CDF 




D0 






all-j 


1+j 


di-1 


1+j 


di-1 


1+j 


di-1 


all-j 


Ixy 


1+j 


di-1 


Lumi (fb-^) 


0.11 


0.11 


0.11 


0.13 


0.13 


0.9 


1.0 


1.0 


0.7 


0.9 


1.0 


Result 


186.0 


176.1 


167.4 


180.1 


168.4 


170.9 


164.5 


171.1 


183.9 


170.5 172.5 


iJES 


0.0 


0.0 


0.0 


0.0 


0.0 


1.4 


0.0 


2.4 


0.0 


0.0 


0.0 


aJES 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.6 


2.0 


bJES 


0.6 


0.6 


0.8 


0.7 


0.7 


0.6 


0.6 


0.4 


0.0 


0.5 


1.8 


cJES 


3.0 


2.7 


2.6 


2.0 


2.0 


0.0 


2.8 


0.0 


0.0 


0.0 


4.3 


dJES 


0.3 


0.7 


0.6 


0.0 


0.0 


0.2 


1.6 


0.0 


0.0 


1.6 


1.9 


rJES 


4.0 


3.4 


2.7 


2.5 


1.1 


0.0 


1.3 


0.0 


0.3 


0.0 


0.0 


Signal 


1.8 


2.6 


2.8 


1.1 


1.8 


1.1 


0.9 


1.3 


1.4 


0.6 


0.7 


BG 


1.7 


1.3 


0.3 


1.0 


1.1 


0.2 


0.7 


1.0 


2.3 


0.3 


0.6 


Fit 


0.6 


0.0 


0.7 


0.6 


1.1 


0.4 


0.9 


0.7 


4.8 


0.4 


0.9 


MC 


0.8 


0.1 


0.6 


0.0 


0.0 


0.2 


0.9 


1.0 


0.7 


0.0 


0.0 


UN/MI 


0.0 


0.0 


0.0 


1.3 


1.3 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Syst. 


5.7 


5.3 


4.9 


3.9 


3.6 


1.9 


3.9 


3.2 


5.6 


2.0 


5.6 


Stat. 


10.0 


5.1 


10.3 


3.6 


12.3 


1.6 


3.9 


2.8 


14.8 


1.8 


5.8 


Total 


11.5 


7.3 


11.4 


5.3 


12.8 


2.5 


5.6 


4.3 


15.8 


2.7 


8.0 



Table 1: Summary of the measurements used to determine the world average Mt. All numbers 
are in GeV/c^. The error categories and their correlations are described in the text. The 
total systematic uncertainty and the total uncertainty are obtained by adding the relevant 
contributions in quadrature. 

several components in order to more accurately accommodate our best estimate of the relevant 
correlations. Each error category is discussed below. 

Statistical: The statistical uncertainty associated with the Mt determination. 

iJES: That part of the JES uncertainty which originates from in-situ calibration procedures 
and is uncorrelated among the measurements. In the combination reported here it cor- 
responds to the statistical uncertainty associated with the JES determination using the 
W qq' invariant mass in the CDF Run-11 1+j and all-h measurements. Residual JES 
uncertainties, which arise from effects not considered in the in-situ calibration, are in- 
cluded in other categories. 
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aJES: That part of the JES uncertainty which originates from differences in detector e/h 
response between 6-jets and hght-quark jets. It is specific to the D0 Run-II measurements 
and is taken to be uncorrelated with the D0 Run-I and CDF measurements. 

bJES: That part of the JES uncertainty which originates from uncertainties specific to the 
modehng of 6-jets and which is correlated across all measurements. For both CDF and 
D0 this includes uncertainties arising from variations in the semi-leptonic branching 
fraction, 6-fragmentation modeling, and differences in the color flow between 6-jets and 
light-quark jets. These were determined from Run-II studies but back-propagated to the 
Run-I measurements, whose rJES uncertainties (see below) were then corrected in order 
to keep the total JES uncertainty constant. 

cJES: That part of the JES uncertainty which originates from modeling uncertainties corre- 
lated across all measurements. Specifically it includes the modeling uncertainties associ- 
ated with light-quark fragmentation and out-of-cone corrections. 

dJES: That part of the JES uncertainty which originates from limitations in the calibra- 
tion data samples used and which is correlated between measurements within the same 
data-taking period (ie. Run I or Run II) but not between experiments. For CDF this 
corresponds to uncertainties associated with the ^^-dependent JES corrections which are 
estimated using di-jet data events. For D0 Run-II this corresponds to uncertainties as- 
sociated with the light-quark response as determined using the W — > qq' invariant mass 
in the 1+j channel and propagated to the di-1 channel. The residual r^-dependent and 
PT-dependent uncertainties for the D0 Run-II measurements are also included here since 
they are constrained using Run-II 7-l-jet data samples. 

rJES: The remaining part of the JES uncertainty which is correlated between all measurements 
of the same experiment independent of data-taking period, but is uncorrelated between 
experiments. This is dominated by uncertainties in the calorimeter response to light- 
quark jets. For CDF this also includes small uncertainties associated with the multiple 
interaction and underlying event corrections. 

Signal: The systematic uncertainty arising from uncertainties in the modeling of the it signal 
which is correlated across all measurements. This inchides uncertainties from variations 
in the ISR, FSR, and PDF descriptions used to generate the tt Monte Carlo samples 
that calibrate each method. It also includes small uncertainties associated with biases 
associated with the identification of 6-jets. 

Background: The systematic uncertainty arising from uncertainties in modeling the domi- 
nant background sources and correlated across all measurements in the same channel. 
These include uncertainties on the background composition and shape. In particular 

uncertainties associated with the modeling of the QCD multi-jet background (all-j and 
l-|-j), uncertainties associated with the modeling of the Drell-Yan background (di-1), and 
uncertainties associated with variations of the fragmentation scale used to model W-|-jets 
background (all channels) are included. 
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Fit: The systematic uncertainty arising from any source specific to a particular fit method, 
including the finite Monte Carlo statistics available to calibrate each method. 



Monte Carlo: The systematic uncertainty associated with variations of the physics model 
used to calibrate the fit methods and correlated across all measurements. For CDF it 
includes variations observed when substituting PYTHIA [291 EOl [31] (Run I and Run II) 
or ISAJET [32] (Run I) for HERWIG [Ml [M] when modehng the tt signal. Similar 
variations are included for the D0 Run-I measurements. The D0 Run-II measurements 
use ALPGEN [35] to model the ti signal and the variations considered are included in 
the Signal category above. 

UN/MI: This is specific to D0 and includes the uncertainty arising from uranium noise in the 
D0 calorimeter and from the multiple interaction corrections to the JES. For D0 Run-I 
these uncertainties were sizable, while for Run-II owing to the shorter integration time 
and in-situ JES determination, these uncertainties are negligible. 

These categories represent the current preliminary understanding of the various sources of 
uncertainty and their correlations. We expect these to evolve as we continue to probe each 
method's sensitivity to the various systematic sources with ever improving precision. Varia- 
tions in the assignment of uncertainties to the error categories, in the back-propagation of the 
bJES uncertainties to Run-I measurements, in the approximations made to symmetrize the 
uncertainties used in the combination, and in the assumed magnitude of the correlations all 
negligibly effect (-C O.lGeV/c^) the combined Mt and total uncertainty. 



4 Correlations 

The following correlations are used when making the combination: 

• The uncertainties in the Statistical, Fit, and iJES categories are taken to be uncorrelated 
among the measurements. 

• The uncertainties in the aJES and dJES categories are taken to be 100% correlated among 
all Run-I and all Run-II measurements on the same experiment, but uncorrelated between 
Run I and Run II and uncorrelated between the experiments. 

• The uncertainties in the rJES and UN/MI categories are taken to be 100% correlated 
among all measurements on the same experiment. 

• The uncertainties in the Background category are taken to be 100% correlated among all 
measurements in the same channel. 
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Run-I published 


Run-II preliminary 


CDF 
1+j di-1 all-j 


D0 
1+j di-1 


CDF 

1+j di-1 all-j Ixy 


D0 
1+j di-1 


CDF-I 1+j 
CDF-I di-1 
CDF-I all-j 


1.00 

0.29 1.00 

0.32 0.19 1.00 








D0-I 1+j 
D0-I di-1 


0.26 0.15 0.14 
0.11 0.08 0.07 


1.00 

0.16 1.00 






CDF-II 1+j 
CDF-II di-1 
CDF-II all-j 
CDF-II Ixy 


0.19 0.13 0.08 
0.36 0.23 0.26 
0.12 0.10 0.10 
0.07 0.03 0.02 


0.14 0.07 
0.24 0.12 
0.08 0.05 
0.05 0.01 


1.00 

0.13 1.00 

0.17 0.10 1.00 

0.05 0.03 0.04 1.00 




D0-II 1+j 
D0-II di-1 


0.12 0.07 0.05 
0.25 0.16 0.17 


0.10 0.04 
0.25 0.11 


0.16 0.06 0.09 0.04 
0.09 0.32 0.05 0.01 


1.00 

0.26 1.00 



Table 2: The resulting matrix of total correlation coefficients used to determined the world 
average top quark mass. 

• The uncertainties in the bJES, cJES, Signal, and Generator categories are taken to be 
100% correlated among all measurements. 

Using the inputs from Tableland the correlations specified here, the resulting matrix of total 
correlation co-efficients is given in Table [2j 

The measurements are combined using a program implementing a numerical minimiza- 
tion as well as the analytic BLUE method [23l[23]. The two methods used are mathematically 
equivalent, and are also equivalent to the method used in an older combination [36], and give 
identical results for the combination. In addition, the BLUE method yields the decomposition 
of the error on the average in terms of the error categories specified for the input measure- 
ments |24j . 



5 Results 



The combined value for the top-quark mass is: 

Mt = 170.9 ± 1.8 GeV/c^ , 
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(1) 





Run-I published 


Run-II preliminary 


CDF 

1+j di-1 all-j 


D0 
1+j di-1 


CDF 

1+j di-1 all-j Ixy 


D0 
1+j di-1 


Pull 

Weight [%] 


+0.73 -0.31 +1.33 
-1.3 -0.4 -0.3 


+ 1.84 -0.20 
+6.1 +0.4 


-0.03 -1.22 +0.05 +0.83 
+39.3 +6.4 +11.0 +0.5 


-0.22 +0.20 
+39.7 -1.9 



Table 3: The pull and weight for each of the inputs used to determine the world average mass 
of the top quark. See Reference [23] for a discussion of negative weights. 



with a of 9.2 for 10 degrees of freedom, which corresponds to a probability of 51% indicating 
good agreement among all the input measurements. The total uncertainty can be sub-divided 
into the contributions from the various error categories as: Statistical (±1.1), total JES (±1.1), 
Signal (±0.9), Background (±0.3), Fit (±0.3), Monte Carlo (±0.2), and UN/MI (±0.1), for a 
total Systematic (±1.5), where all numbers are in units of GeV/c^. The pull and weight for each 
of the inputs are listed in Table [31 The input measurements and the resulting world average 
mass of the top quark are summarized in Figure [H 

The weights of many of the Run-I measurements are negative. In general, this situation 
can occur if the correlation between two measurements is larger than the ratio of their total 
uncertainties. This is indeed the case here. In these instances the less precise measurement 
will usually acquire a negative weight. While a weight of zero means that a particular input 
is effectively ignored in the combination, a negative weight means that it affects the resulting 
central value and helps reduce the total uncertainty. See reference [23] for further discussion of 
negative weights. 

Although the from the combination of all measurements indicates that there is good 
agreement among them, and no input has an anomalously large pull, it is still interesting to 
also fit for the top-quark mass in the all-j, 1+j, and di-1 channels separately. We use the same 
methodology, inputs, error categories, and correlations as described above, but fit for the three 
physical observables, Mf^~'\ mI^\ and Mf^~^. The results of this combination are shown in 
Table [Hand have of 5.8 for 8 degrees of freedom, which corresponds to a probability of 60%. 
These results differ from a naive combination, where only the measurements in a given channel 
contribute to the determination in that channel, since the combination here fully accounts 
for all correlations, including those which cross-correlate the different channels. Using the 
results of Tabled] we calculate the chi-squared consistency between any two channels, including 
all correlations, as x^{.dil — Ij) = 3.2, x^(^J ~ allj) = 0.1, and x^i.'^^^j ~ dil) = 2.4. These 
correspond to chi-squared probabilities of 7%, 75%, and 12%, respectively, and indicate that 
the determinations of from the three channels are consistent with one another. 
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Parameter 


Value (GeV/c^) 


Correlations 




172.2 ±4.1 
171.2 ± 1.9 
163.5 ± 4.5 


1.00 

0.21 1.00 

0.15 0.30 1.00 



Table 4: Summary of the combination of the nine measurements by CDF and D0 in terms of 
three physical quantities, the mass of the top quark in the all-jets, lepton+jets, and di-lepton 
channel. 

6 Summary 

A preliminary combination of measurements of the mass of the top quark from the Tevatron 
experiments CDF and D0 is presented. The combination includes five published Run-I mea- 
surements and two published plus four preliminary Run-II measurements. Taking into account 
the statistical and systematic uncertainties and their correlations, the preliminary world-average 
result is: = 170.9 ± 1.8 GeV/c^, where the total uncertainty is obtained assuming Gaussian 
systematic uncertainties and adding them plus the statistical uncertainty in quadrature. 

The mass of the top quark is now known with a relative precision of 1.1%, limited by 
the systematic uncertainties, which are dominated by the jet energy scale uncertainty. This 

systematic is expected to improve as larger data sets are collected since new analysis techniques 
constrain the jet energy scale using in-situ W — * qq' decays. It can be reasonably expected 
that with the full Run-II data set the top-quark mass will be known to much better than 
1%. To reach this level of precision further work is required to determine more accurately the 
various correlations present, and to understand more precisely the 6-jet modeling. Signal, and 
Background uncertainties which may limit the sensitivity at larger data sets. Limitations of 
the Monte Carlo generators used to calibrate each fit method may also become important as 
the precision reaches the ~ 1 GeV/ level and will warrant further study in the near future. 
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Mass of the Top Quark (*Preliminary) 
Measurement 



Tevatron Run 



M^^p [GeV/cl 




X /dot = 9.2/10 



170.9± 1.8 



150 170 

M^^p [GeV/c'] 



190 



Figure 1: A summary of the input measurements and resulting world average mass of the top 
quark. 
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